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ABSTRACT 

A moderately spinning neutron star acquires an oblate shape and a spacetime with a significant 
quadrupole moment. These two properties affect its apparent surface area for an observer at infinity, 
as well as the light curve arising from a hot spot on its surface. In this paper, we develop a ray-tracing 
algorithm to calculate the apparent surface areas of moderately spinning neutron stars making use 
of the Hartle-Thorne metric. This analytic metric allows us to calculate various observables of the 
neutron star in a way that depends only on its macroscopic properties and not on the details of 
its equation of state. We use this algorithm to calculate the changes in the apparent surface area, 
which could play a role in measurements of neutron-star radii and, therefore, in constraining their 
equation of state. We show that whether a spinning neutron star appears larger or smaller than its 
non-rotating counterpart depends primarily on its equatorial radius. For neutron stars with radii 
~ 10 km, the corrections to the Schwarzschild spacetime cause the apparent surface area to increase 
with spin frequency. In contrast, for neutron stars with radii ~ 15 km, the oblateness of the star 
dominates the spacetime corrections and causes the apparent surface area to decrease with increasing 
spin frequency. In all cases, the change in the apparent geometric surface area for the range of observed 
spin frequencies is < 5% and hence only a small source of error in the measurement of neutron-star 
radii. 

Subject headings: gravitation — relativistic processes — stars: neutron — stars: rotation 



1. INTRODUCTION 

Spectroscopic measurements of thermally emitting 
neutron-star surface areas have been widely used in at- 
tempts to infer the stellar radii and constrain the equa- 
tion of state of ultradense matter (e.g., van Paradijs 1978, 
1979; Rutledge et al. 2001; Majczyna & Madej 2005; 
Ozel 2006). This approach has recently flourished with 
the use of high quality spectra obtained with current X- 
ray satellites during thermonuclear flashes (Ozel et al. 
2009; Giiver et al. 2010a, 2010b) and in the quiescent 
phases of accreting neutron stars (Heinke et al. 2006; 
Webb & Barret 2007; Guillot et al. 2011). In some cases, 
the measured apparent surface areas possess formal un- 
certainties as small as a few percent (see, e.g., Giiver 
et al. 2011) making it possible to distinguish between 
the predictions of different equations of state (Ozel et al. 
2010; Steiner et al. 2010). 

Converting a measurement of the apparent surface area 
of a neutron star to a radius requires correcting for a 
number of general-relativistic effects. As the photons 
propagate in the strong gravitational fields of neutron 
stars, their energies are redshifted while their trajecto- 
ries experience strong lensing. For a non-spinning neu- 
tron star, both effects depend only on the tt-component 
of its exterior metric and can be analytically corrected 
for in any metric theory of gravity (Psaltis 2008). The 
situation, however, becomes increasingly complicated as 
the angular velocity of the neutron star increases toward 
the point of breakup. 

The qualitative character of the general-relativistic ef- 
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of the spin frequency of the neutron star, /g, and the 
characteristic spin frequency fo = \/ GM/R^ /2tt (Hartle 
& Thorne 1968), where G is the gravitational constant 
and M and R are the mass and radius of the neutron 
star, respectively. To zeroth order in /s//o, the neutron 
star is spherically symmetric and its external spacetime 
is described by the Schwarzschild metric, which depends 
only on the mass of the star. To first order in fs/ fo, 
the neutron star remains spherically symmetric and its 
external spacetime is that of the Kerr metric. In this 
case, the amount of gravitational lensing depends also 
on the spin angular momentum of the star (because of 
the effects of frame dragging), which in turn depends 
on the density profile inside the star and hence on the 
equation of state. To second order in fs/ fo, the neu- 
tron star becomes oblate and its external spacetime is 
described by the Hartle-Thorne metric, which depends 
on the mass of the neutron star, on its spin angular mo- 
mentum, and on its quadrupole mass moment. Finally, 
to even higher orders in fs/ fo, the external spacetime 
of the neutron star depends on multipole moments that 
are of increasing order. These spacetimes can be calcu- 
lated numerically (see, e.g., Cook et al. 1994; Stergioulas 
& Friedman 1995). In principle the external spacetime 
of a rapidly spinning neutron star can be accurately de- 
scribed by the analytic solution of Manko et al. (2000a, 
2000b; see Berti & Stergioulas 2004; Berti et al. 2005). 
However, the Manko et al. (2000a, 2000b) metric does 
not reduce to the Schwarzschild or Kerr solutions when 
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the angular velocity of the star is reduced toward zero 
(Berti &: Stergioulas 2004), making it impractical for ap- 
plications. 

Various effects of increasing the angular velocity of a 
neutron star have been studied in the context of predict- 
ing the orbits of test particles (e.g., Shibata & Sasaki 
1998; Abramowicz et al. 2003; Berti & Stergioulas 2004; 
Berti et al. 2005), the lightcurves that arise when the sur- 
face emission on a spinning neutron star is not uniform 
(Miller & Lamb 1998; Braje et al. 2000; Muno et al. 2003; 
Poutanen & Gicrhnski 2003; Cadeau et al. 2005, 2007; 
Morsink et al. 2007), as well as the rotational broaden- 
ing of atomic lines that originate on the stellar surfaces 
(Ozel & Psaltis 2003; Bhattacharya et al. 2006; Chang et 
al. 2006). There are two conclusions that emerged from 
these studies and are relevant for our work: (1) that the 
Hartle-Thorne metric provides an approximation to the 
external spacetime of a rotating neutron star that is ad- 
equate for most astrophysical applications (Berti et al. 
2005) and (2) that taking into account the oblateness of 
the stellar surface is at least as important as considering 
the effects of frame dragging around a spinning neutron 
star (Morsink et al. 2007). 

In this article, we calculate numerically the appar- 
ent geometric surface area of a rotating neutron star 
using a variant of the Hartle-Thorne metric developed 
by Glampedakis & Babak (2006). We use our new, 
fast ray-tracing algorithm that takes into account the 
oblateness of the neutron star as well as deviations 
from the Schwarzschild metric for its external spacetime 
that are formally correct up to the contributions of the 
quadrupole mass moment. With this analytic metric we 
can calculate the apparent surface area of a neutron star 
in a way that depends only on global properties of the 
star and not on the details of an assumed equation of 
state. Indeed, as we will show, the apparent surface area 
depends on six parameters: the equatorial and polar radii 
of the neutron star, its mass, spin angular momentum, 
and quadrupole mass moment, as well as the observer's 
inclination. Using the approximate relation between the 
polar and equatorial radii of spinning neutron stars ob- 
tained by Morsink et al. (2007) reduces the number of 
parameters to five. 

We find that the apparent surface area of the neutron 
star changes only marginally over the range of spin fre- 
quencies of known sources. The corrections to the ob- 
served surface area only become significant at rotational 
frequencies greater than 1000 Hz. 

2. THE SPACETIME OF A RAPIDLY SPINNING NEUTRON 

STAR 

We describe the external spacetime of a rapidly spin- 
ning neutron star using the variant of the Hartle-Thorne 
metric developed by Glampedakis & Babak (2006). This 
is based on the Kerr metric, but with the quadrupole mo- 
ment allowed to deviate from its Kerr value. It therefore 
describes the spacetime of a rapidly spinning neutron star 
up to second order in /s//o- We opted to use this metric 
in order to contrast our results with those obtained when 
the external spacetime of a rapidly spinning neutron star 
is approximated by the Kerr metric. 

In Boyer-Lindquist coordinates, we write the metric as 
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The quadrupole correction is given, in contravariant 
form, by 
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with the functions J^i,2('') shown explicitly in Ap- 
pendix A of Glampedakis & Babak (2006) . In these equa- 
tions and hereafter, we set G = c = 1. Note that, because 
it is of even order, the mass quadrupole affects only the 
diagonal components of the Kerr metric and not the tcj)- 
component that measures the amount of frame dragging. 
In this approach, the apparent surface area of a spinning 
neutron star depends on at least three macroscopic prop- 
erties of the spacetime: the mass M, the specific angular 
momentum per unit mass, a, and the mass quadrupole 
moment q. 
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Fig. 1. — Mass quadrupole moment of the spacetime as a function 
of the spin frequency of the neutron star for two different radii. For 
each configuration, we show the quadrupole moment for a star of 
mass 1.8 Mq, with a radius of 10 km and 15 km, and the remaining 
parameters chosen to correspond to models with EOS EPS and 
EOS L, respectively (Cook et al. 1994). 



Images of Spinning Neutron Stars 



3 



For a compact object with mass M and specific angular 
momentum a, we write the mass quadrupole moment of 
the spacetime as 



q = -a''{l + ri) , 



(7) 



so that, when rj ~ 0, this reduces to the quadrupole mo- 
ment of the Kerr metric. (Note that, in the formalism of 
Glampedakis & Babak 2006, e = rja^.) This mass mo- 
ment of the spacetime for a rapidly spinning neutron star 
depends on its density profile and hence on the underly- 
ing equation of state. 

Laarakkers & Poisson (1999) calculated neutron-star 
quadrupole moments for a wide range of equations of 
state, using a numerical algorithm that solves the Ein- 
stein field equations with no approximations. They found 
that Equation ([7|) remains valid even for neutron stars 
that are spinning near their breakup points. They cal- 
culated 77 ^ 1 — 6, depending on the stellar mass and 
radius. Berti & Stergioulas (2004) further computed the 
quadrupole moments of stars in the Hartle-Thorne ap- 
proximation and compared them to those of the numeri- 
cal spacetimes. They concluded that the Hartle-Thorne 
approximation is adequate for all astrophysical applica- 
tions. This approximation is expected to be valid for 
observed neutron stars because the fastest known X-ray 
burster is spinning at 620 Hz (4U 1608-52; Galloway et 
al. 2008), and the fastest known pulsar is spinning at 
716 Hz (Hessels et al. 2006), which are expected to be 
significantly smaller than their breakup frequencies. Fig- 
ure [1] shows the dependence of the quadrupole moment q 
on the spin frequency for a 1.8 Mq star with the appro- 
priate value of the parameter rj taken from Laarakkers & 
Poisson (1999) for two different equations of state. 

In order to calculate the apparent surface area of a 
neutron star using the above spacetime, we need to al- 
low for the stellar surface to be non-spherical. This is re- 
quired by the fact that the deviations of both the external 
spacetime from the Kerr metric as well as of the shape 
of the stellar surface from spherical symmetry are of sec- 
ond order in fs/ fa- Moreover, Morsink et al. (2007) have 
shown that, for purely geometric reasons, the changes in 
the predicted lightcurves of spinning neutron stars when 
the non-spherical shapes of their surfaces are taken into 
account are at least as important as the effects of frame 
dragging that are formally only of first order in /s//o. 

Morsink et al. (2007) fit a large number of neutron-star 
shapes calculated for a wide range of equations of state. 
They show that an equation of the form 
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Fig. 2. — Ratio of tiie polar to tiie equatorial radius of a neutron 
star as a function of its spin frequency, calculated using the approx- 
imate relation of Morsink et al. (2007). The two pairs of curves 
correspond to two different stars with the same mass of 1.8 Mq 
and different radii. In each case, the solid curve is the result of 
the complete fitting formula whereas the dashed curve is the 
result when only the terms up to the quadrupole order are used. 
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Excluding those equations of state that describe self- 
bound strange stars, Morsink et al. (2007) obtain for 
the coefficients in Equation (|5]) the following empirical 
relations: 



ao = -0.18e 
a2 = -0.39e 
a4 = 0.04e-0.15Ce-f 0.07e 



0.23Ce- O.OSe^ 
0.29Ce-H0.13e2 



(11) 



Figure [5] shows the ratio of the polar to the equatorial 
radius of a neutron star, as calculated with the above 
fitting formula. 

In the Hartle-Thorne approximation, only terms up 
to second order are considered, i.e., the above series is 
terminated at A'^ = 1. In this case, we can rewrite the 
shape of the stellar surface in terms only of the equatorial 
and polar radii as 
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with i?cq the equatorial radius of the star and P2«(cos6') This introduces the two additional parameters, i?p and 

i?cq, which will determine the apparent surface area of 
the neutron star. In Figure [21 we show the deformation 
of neutron stars for both the full fitting formula (solid 
lines) and the approximation up to quadrupole order 
(dashed lines). Over the entire region of interest, the 
agreement between the approximation and the exact for- 
mula is within 5%. Therefore disregarding higher order 
mass moment terms when calculating the shape of the 



l-cq 

the Legendre polynomial of order 2n, accurately de- 
scribes the shape of the neutron-star surface for even 
the fastest spinning stars if the series is terminated at 
N — 2. For all equations of state the coefficients of the 
expansion depend only on two parameters. 
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Configuration 1 



Configuration 2 




Fig. 3. — Apparent image of neutron stars at infinity calculated 
for different configurations. Two stars witii Rcq = 15 km (left pan- 
els) and Rcq = 10 km (right panels) are shown, one non-spinning 
(upper panels) and one with a spin frequency of 1400 Hz (lower 
panels). In both cases, the observer is in the plane of the equa- 
tor {00 = ^), and the star rotates counterclockwise, the left side 
approaching the observer and the right side receding. 

neutron star is justified. 

The last free parameter in our simulation is the incli- 
nation angle of the observer. Since there is no spherical 
symmetry, the apparent surface area of the neutron star 
will depend on the angle at which it is viewed. We define 
this angle, 9o, as the angle between the rotation axis of 
the star and the line of sight to an observer at infinity. 

3. A NUMERICAL ALGORITHM FOR RAY TRACING 

We follow the method outlined in Psaltis & Johannsen 
(2011) to develop a ray-tracing algorithm that calculates 
the trajectories of photons in the neutron-star spacetime. 
The algorithm integrates the geodesic equations of the 
Hartle-Thorne metric using a Runge-Kutta solver. A 
number of rays originate perpendicular to an image plane 
at a large distance from the neutron star and are traced 
backward until they either intersect the surface of the 
star or exceed the distance to the star by 10%, at which 
point we assume that they have bypassed the star en- 
tirely. 

Because the surface of the neutron star deviates from 
spherical symmetry at high spin frequencies and photons 
may intersect the surface at oblique angles, special care 
must be taken when evaluating the coordinates at which 
the photon geodesies intersect the stellar surface. We 
use a simple procedure to determine the spherical coor- 
dinates of intersection: once a step in the Runge-Kutta 
algorithm crosses over the surface of the star, we perform 
a linear interpolation between the last two points on the 
geodesic. We then employ a bisection algorithm to find 
the point of intersection between that interpolation and 
the functional form of the stellar surface. 

Following Cadeau et al. (2007), we calculate the beam- 



ing angle for a photon emitted at the surface of the neu- 
tron star as 



COSde = 



k"- - k^R'{9) 
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where k"^ and k^ are the components of the photon's mo- 
mentum at the point of intersection with the stellar sur- 
face, R{9) is the functional form of the surface of the 
neutron star (Equation ([5])), and z is the total redshift. 

4. GEOMETRIC SURFACE AREAS OF NEUTRON STARS 

One way of inferring the apparent surface area of the 
neutron star, in principle, is by direct angular measure- 
ment of the stellar image. If the distance of the star 
is known, it is straightforward to calculate its physical 
radius given its angular dimension. Another method to 
infer the surface area is by means of its thermal spec- 
trum. By measuring the temperature of the star and 
its total flux, one can again derive its surface area. For 
slowly spinning neutron stars, these two measurements 
agree. 

In the absence of rotation, the spacetime around the 
star is described by the Schwarzschild metric. An ob- 
server at infinity will observe the star to be enlarged due 
to gravitational self-lensing. The angular apparent sur- 
face area of the star is then given by 



4^ 
A ^ 



2GM 
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where A = ttR^. The spectroscopic area, on the other 
hand, is defined as 



4'S' 

^Sch 



(15) 



where Foo is the thermal fiux and Tcff.oo is its effective 
temperature. Gravitational redshift reduces both the ob- 
served flux and the effective temperature as 

2GM 



Rc^ 



2GM 
1 ^t^Tns, 



(16) 



where -Fns ^-nd Tns are the flux and temperature mea- 
sured at the surface of the neutron star. Combining 
Equations (fT5)) and ([T5)) we find that the apparent spec- 
troscopic surface area at infinity is given by 



"^Sch 



4(1 



2GM 



A i?c2 ' ■ 

Clearly, for a slowly spinning neutron star 

^Sch _ A 
^Sch 



(17) 



(18) 



The factor of four arises from the difference between the 
total surface area (47ri?^) measured spectroscopically and 
the projected surface area {ttR^) measured geometrically. 
For a moderately spinning neutron star, the lack of spher- 
ical symmetry in the metric, the oblateness of the star. 
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Fig. 4. — Effect of frame dragging, quadrupole moment, and oblateness on the image of a neutron star. The left panel shows the image 
of a star in a Kerr metric with q = 0, where the only effect on the appearance of the star is the displacement of the image due to frame 
dragging. The center panel shows the star in a Kerr metric with an extra quadrupole moment as described in Equation (0, which causes 
the image to stretch horizontally. The right panel includes both an extra quadrupole and an oblate stellar surface, causing the poles to 
squeeze together. The mass, radius, and moment of inertia for each star are chosen according to Configuration 2, as is the quadrupole 
moment for all but the Kerr case. In each case the star is rotating at 1000 Hz. The red dashed circle in each panel indicates the location 
and size of the image in the Schwarzschild metric. 



TABLE 1 

Neutron- Star Parameters 



Configuration 


M 




/ 


V 




(Me) 


(km) 


Clod's gcm^) 




1 


1.8 


10 


1.4 


1.2 


2 


1.8 


15 


3.0 


3.9 



Parameters of interest for the two neutron star configurations used 
in this paper: the mass A/, the equatorial radius -Rcq, the moment 
of inertia /, and the deviation from the Kerr quadrupole rj. 



and the second order terms in the Doppler shift will in- 
troduce different corrections to the geometric and spec- 
troscopic areas as measured at infinity. As a result, these 
two measurements will not agree in general. 

In practice, only the spectroscopic area can be mea- 
sured, as neutron stars are too small to allow for a di- 
rect measurement of their angular sizes. However, the 
change in apparent angular size of a neutron star plays a 
role even in the measurement of its spectroscopic surface 
area. In the following, we use our ray-tracing algorithm 
to investigate the effect of rapid rotation on the geometric 
surface areas of neutron stars. We will discuss the addi- 
tional effects of position-dependent redshift and Doppler 
shifts in this metric in a forthcoming paper. 

We calculate the apparent geometric surface area of 
a neutron star by counting directly the number of rays 
originating at the image plane that intersect the surface 
of the star. The area on the image plane covered by the 
intersecting rays corresponds to the angular area of the 
image of the neutron star. Moreover, the coordinates on 
the neutron star at which each ray intersects the surface 
allow us to construct a contour plot of the apparent image 



of the star to an observer at infinity (see Figure |3]). 

For our models, we used two different configurations 
of neutron-star parameters, hereafter labeled "Configura- 
tion 1" and "Configuration 2" . In each case, we chose the 
mass, equatorial radius, moment of inertia (from Cook 
et al. 1994), and quadrupole deviation r] (according to 
Laarakkers & Poisson 1999) corresponding to a proposed 
equation of state in the non-spinning case. In order to 
demonstrate the range of the change in surface area, we 
chose a relatively soft equation of state (EOS FPS) for 
Configuration 1, and a relatively stiff equation of state 
(EOS L) for Configuration 2. We then held these pa- 
rameters constant and varied only the spin frequency of 
the neutron star. Therefore, as the spin of each star in- 
creases, its parameters will increasingly deviate from the 
values predicted by the initial equations of state. While 
the resulting models do not reflect the predictions of any 
particular equation of state, they allow us to separate 
the various physical effects and estimate the magnitude 
of the change in surface area across the parameter space. 
The parameters of both configurations are summarized 
in Table [11 

Figure [3] illustrates the appearance of the two neutron 
stars in both the Schwarzschild approximation and at a 
high spin frequency. Three effects influence the appar- 
ent image in the rapidly spinning case: frame dragging, 
the non-zero quadrupole moment, and the oblateness of 
the stellar surface. The change in the apparent image 
of the star due to each of these effects is illustrated in 
Figured! Frame dragging causes the image of the neu- 
tron star to be displaced horizontally while preserving 
its circular appearance. The negative quadrupole mo- 
ment causes the image to stretch horizontally about the 
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Fig. 5. — Ratio of the apparent surface area of a rotating neutron star to the area of a static neutron star as a function of spin frequency. 
Two stars are shown, corresponding to Configuration 1 (left panel) and Configuration 2 (right panel) as described in the text. In both plots, 
the four lines show different approximations of the metric outside the star: the Kerr metric, the Kerr metric plus a quadrupole deviation, 
the Kerr metric with an oblate stellar surface, and the Kerr metric with both a quadrupole deviation and oblateness. The dashed line is 
at A/Asch = 1- In this plot, the inclination of the star is set to 90°. Note that the plot is extended past the mass shed limit in some 
configurations in order to show the general trend of the relationship. 




^•^ 500 1000 1500 500 1000 1500 
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Fig. 6. — As in Figure[5]but showing the apparent radius for three different inclinations, taking into account both the quadrupole moment 
and the oblateness of the star. The left panel corresponds to Configuration 1, while the right panel corresponds to Configuration 2. Again, 
the X-axis is extended past the mass-shed limit in some cases. 



axis of rotation. Finally, the oblateness of the stellar sur- 
face squeezes the poles together, further emphasizing the 
elongated appearance of the star. The combination of the 
displacement due to frame dragging and the stretching 
due to the quadrupole moment give rise to the asymmet- 
ric shape of the stellar image. A similar effect has been 
reported for images of accretion flows around black holes 
(Johannsen & Psaltis 2010). 

Figure O shows the contributions of each of these ef- 
fects to the geometric surface area of the neutron star 
as a function of spin frequency. Note that the Hartle- 
Thorne approximation to the metric formally becomes 
inaccurate at the highest spin frequencies shown in this 



figure. The dominant effect on the image of a neutron 
star with a Kerr metric is due to frame dragging. How- 
ever, as discussed above, frame dragging merely displaces 
the image of the star without altering its area. Therefore, 
the apparent area under the Kerr metric stays nearly con- 
stant (the small increase at high spin frequencies is due 
to the non-zero Kerr quadrupole). Including the extra 
quadrupole moment described by Equation (jT)) increases 
the surface area more significantly at faster spins. The 
oblateness of the stellar surface tends to counteract these 
effects and decrease the apparent surface area. Whether 
the surface area increases or decreases with higher spin 
frequency depends on which of these corrections domi- 
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Fig. 7. — Relative contributions of geometry and gravitational 
light bending to the total surface area as a function of inclination. 
The green line, marked "Geometry," shows the apparent surface 
area normalized to ttR^^ in a fiat spacetime (see the Appendix). 
The red line, labeled "Total," shows the area calculated by ray 
tracing in the Hartle-Thorne metric. The ratio of these two ar- 
eas shows the effect of gravitational light-bending as a function of 
inclination, plotted in the blue line labeled "Gravity" . 

nate. 

Comparing the left and right panels of Figure [51 it is 
clear that the relative importance of the oblate shape of 
the star grows as its equatorial radius increases. Fig- 
ure [5] shows that the radius of a star described by Con- 
figuration 2 falls within the regime where the oblateness 
dominates the correction to the apparent surface area. 
In this case, the observed surface area of the star pri- 
marily decreases with increasing spin frequency, and the 
quadrupole moment only adds a small correction at the 
highest spin frequencies. At a radius of 10 km, however, 
as in Configuration 1, the effect of the quadrupole mo- 
ment outweighs the oblateness for spin frequencies up 
to 1200 Hz, resulting in an overall increase of apparent 
surface area. 

Since the presence of a spin axis breaks the symmetry 
of the spacetime, the image of the neutron star also de- 
pends on the inclination at which it is viewed. Figure |6] 
shows this dependence for the two neutron stars with dif- 
ferent radii. Again, the direction of the trend depends 
on which of the effects is dominant. For the neutron star 
with a larger radius (left panel), increasing the spin in- 
variably leads to a smaller observed surface area and a 
weak dependence on inclination. For the neutron star 
with a smaller radius, however, the effects due to the 
oblateness and the quadrupole moment are comparable 
in magnitude (right panel). In this case, the direction of 
the effect depends more strongly on inclination: the star 
appears to grow in area with increasing spin frequency 
at large inclination, while it shrinks when viewed near 
the poles. 

The change in surface area with inclination is due to 
the non-spherical shape of the neutron star and the 6- 
dependence of the metric. We investigate the behavior 
of these two contributions further in Figure [71 where we 



show the contributions of geometry and light-bending to 
this change in area for a fixed spin frequency for the 
Configuration 2 neutron star. 

The curve labeled "Geometry" shows the dependence 
of the apparent surface area of the neutron star on in- 
clination in the absence of any relativistic effects (see 
the Appendix for the details of the calculation). As ex- 
pected, changing the inclination from pole-on to edge-on 
decreases the surface area of the image. On the other 
hand, the curve labeled "Total" shows the change of the 
surface area with inclination when relativistic effects are 
taken into account. We also plot the ratio of the two de- 
pendences as the curve labeled "Gravity." As the poles 
of the oblate spheroid lie deeper in the gravitational well 
than the equator, the effect of light-bending is stronger 
when the star is seen at a higher inclination. These two 
effects cancel to within 10%, leaving the apparent sur- 
face area of the star practically constant over changing 
inclination. 

5. CONCLUSIONS 

A spinning neutron star introduces several effects in 
surface area measurements beyond those of the station- 
ary Schwarzschild metric. Rapid rotation causes the neu- 
tron star to acquire a quadrupole moment and an oblate 
surface. Both of these effects distort the image of the star 
as seen by distant observers, introducing a correction to 
the calculation of neutron star radii from a measurement 
of their surface areas. 

The relative contributions of the quadrupole moment 
compared to the oblateness depend strongly on the ra- 
dius of the star. At a radius of 15 km, the quadrupole 
contribution is negligible when compared to the change 
in apparent area caused by the deformation of the stellar 
surface at all spin frequencies. At a radius of 10 km, 
which are consistent with current observations (Ozel 
et al. 2010; Steiner et al. 2010), the self-lensing domi- 
nates over the oblateness for spin frequencies less than 
- 1000 Hz. 

While the apparent surface area does depend on 
the oblateness and quadrupole moment caused by the 
neutron-star spin, these corrections are small for the spin 
frequencies that have been observed to date. Only at fre- 
quencies above ~ 800 Hz does the apparent surface area 
change significantly. For the observed range of spin fre- 
quencies, the change in apparent geometric surface area 
is less than 10%, for both configurations considered here 
and at all inclinations. 

Studies of the spectroscopic determination of surface 
areas for thermally emitting neutron stars have so far re- 
lied on the Schwarzschild approximation, which neglects 
the correction due to the changes in spin and oblateness 
of the neutron star due to the stellar rotation (e.g.. Rut- 
ledge et al. 2001; Ozel 2006; Heinke et al. 2006; Webb & 
Barret 2007; Ozel et al. 2009; Guver et al. 2010a, 2010b). 
Observational uncertainties inherent in these measure- 
ments are currently at the ^ 10% level and hence larger 
than the effects presented here. We will discuss the addi- 
tional corrections in the spectroscopic measurements of 
spinning neutron stars in a forthcoming paper. As the 
uncertainties in the data shrink with future observations, 
it will be necessary to take into account the corrections 
due to neutron star spins. 
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APPENDIX 

In this appendix we calculate the Newtonian surface area of the neutron star as a function of the inclination angle. 
Qualitatively, we are calculating the area of the shadow cast by the neutron star onto the image plane in the presence 
of light rays normal to the image plane. The outline of this figure is defined by the points on the neutron star at which 
a normal to the image plane is tangent to its surface. 

To find this locus of points, we first write a parametric expression for a vector from the origin to a point on the 
surface of the neutron star with coordinates {0,(j)) in the usual form: 

R(6l, 0) = {R{e) sm{e) cos(0), R{0) sm{e) sin((^), R{e) cos{e)}, (Al) 

with R{d) defined as in Equation ((S)). We then define a normal vector to the surface at (6',(/)) as 

and a unit vector normal to the image plane at an angle to the rotational axis of the neutron star as V(0o) = 
{0, sin(0o); cos(0o)}- As 6*0 varies from to ^, the inclination of the image plane ranges from pole-on to edge-on. 

The boundary of the projection onto the image plane is defined by those points on the stellar surface which satisfy 
the equation 

N(0,0)-V(0o) = O. (A3) 
This equation is linear in sin((/)) and has the solution 

. cos(6l)cot(e'o)i?(e') +cot(6'o)sin(6l)i?'(6i) 

= R{e)sMe)-cosie)R'ie) ' ^^'^ 

This equation is no longer valid in the limiting cases — ^ and 6 = 0. In the former case, the outline of the neutron 
star is described by Equation ([8]). In the latter case, when the observer lies along the rotational axis of the star, the 
outline will of course appear as a circle with radius i?(-|). It should also be noted that not all values of colatitude 
9 will correspond to a physical solution for the azimuth (f> — the projected circumference of the image will lie within 
a range of 9 around the equator. We denote by 9+ and 6^ the limits of colatitude that contribute to the image at 
infinity. 

In order to project the resulting outline of the neutron star onto the image plane, we simply rotate the coordinate 
system by an angle —6*0 so the x-y plane is parallel to the image plane. We then set the z-component of all the points 
on the outline to zero to find a two-dimensional parametric form {Cx{9; 9o), Cy{9; 9o)} of the outline of the projected 
image. Finally, we numerically integrate this curve to find the projected surface area: 

A = 2rC,^d9. (A5) 
. f> d9 
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